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the multi-cut two-matrix models, based on an extension of the prescription de- 
veloped by Daul, Kazakov and Rostov. After identifying possible critical points 
and potentials in the multi-cut matrix models, we calculate the macroscopic loop 
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the matrix Lax operators, a sequence of new solutions for the amplitudes in the Z/% 
symmetric fc-cut two-matrix models are obtained, which are realized by the Jacobi 
polynomials. 
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1 Introduction and summary 



Non-critical string theory [1-3] provides a simple and tractable toy model to understand 
various aspects of string theory. One of the fascinating formulation of the theory is given 
by the description with the matrix model. This realizes a non-perturbative definition of 
string theory [4-6] and provides an explicit example of gauge/string duality [7,8]. 

Recently, our knowledge about the duality in non-critical string theory has extended to 
various kind of string theories. The first non-trivial extension was to the correspondence 
between the two-cut matrix models [9-14] and type superstring theory [15-17]; later, 
the multi-cut matrix models [18] were proposed as a description of fractional superstring 
theory [19]. Despite of many qualitative discussions, little is known quantitatively about 
the nature of multi-cut matrix models. Therefore, in this paper we try to give explicit 
answers to the following questions: how do critical points and potentials look like and 
what is the corresponding geometry appearing at the vicinity of critical points? 

To address these issues, macroscopic loop amplitudes [20-27] have played important 
roles. The one-point function of macroscopic loop operator, which is also known as 
resolvent, is defined as 

Geometrically, the macroscopic loop amplitudes describe the probability amplitudes for 
creating holes on the two-dimensional surface as induced from the Feynman graph of 
matrix integral. Analytically, we can view the macroscopic loop amplitudes as a gener- 
ating function for correlation functions of local operators (in X) in the matrix models. 
Most importantly, it enables us to probe eigenvalue distribution of a matrix X [20] and 
spacetime geometry in the weak coupling region [28,29]. Thus, the macroscopic loop 
amplitude is of fundamental importance in the theory of matrix models. In order to 
calculate the macroscopic loop amplitudes, one powerful way is to use the method of or- 
thogonal polynomials [27]. Through a useful identity (see Appendix , we can express 
the resolvent Q and the spectral parameter x as a pair of conjugate variables acting on the 
space of orthogonal polynomials, whose commutations relation gives rise to the Douglas 
(string) equation. In terms of the differential operators, the string equation turns into 
a differential equation on the macroscopic loop amplitudes and we can solve it without 
directly evaluating the matrix integral (11.11) . 

Right at the critical points, the resolvent has specific critical behavior. Conversely, 
with some proper choice of the critical resolvent, one can obtain the critical potentials. In 
the two-matrix-model cases, this approach is especially useful for finding critical points 
[30,31] and has been explicitly formulated for general (p, q) critical points [27]. Along this 
line, we explicitly construct the critical points and potentials in the multi-cut two-matrix 
models. Our results provide further insight into the understanding of the relationship 
[14] between the existence of critical points and the the choice of hermiticity of critical 
potentials. 

Macroscopic loop amplitudes at the vicinity of critical points (i.e. with perturbation 
of the cosmological constant) have been investigated in [22-24,26,27,32]. The seminal 
formula for general (p, q) critical points with cosmological constant /i was first found by 
Kostov [22] and can be expressed as 

x = y/jiT p (z), Q = ^*T q (z), (1.2) 
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with Tp(coshr) = coshpr, the Chebyshev polynomials of the first kind, and z = coshr. In 
the case of (p, q) minimal superstring theory^] the corresponding formula was also found 
on the Liouville side [34]. The non-trivial solution^] can be expressed^] as 

x = ^Up(z) V^l, Q = fl 7 * Ug(z) v^^l, (1.3) 

with [/p(coshr) = sinhpr/ sinh r, the Chebyshev polynomials of the second kind, and 
z = coshr. This general formula has been utilized in various studies of (p,q) (or (p,q)) 
critical points [32-39], even in the D-instanton effects of the string theory [28,40-43]. 
Moreover, in comparison with worldsheet CFT, e.g. with the FZZT-brane or ZZ-brane 
amplitudes in Liouville theory [44,45], it is essentially the above formula which is realized 
in the calculation of (p,q) (or (p,q)) minimal (super)string theory [34]. 

One of the significances of this formula is that it provides us further checks of the 
string duality (e.g. matching of correlators) beyond matching of operator contents and 
critical exponents. Since this comparison is necessary especially in the correspondence 
between the multi-cut matrix models and fractional superstring theory [19], discovery of 
the general formula for the multi-cut two-matrix models is one of the most important 
stage in understanding the conjectured duality as well as the system itself. 

In this paper, the generalization of these formulae are explored especially in the 
symmetric potential w(x,y), 

w(u n x, u- n y) = w{x, y), u = e 2m/k , (n e Z), (1.4) 

of the general fc-cut two-matrix models (k = 3,4, • • •). Identifying a proper ansatz for 
the non-polynomial part of the amplitudes (generalizing the square-root of (jl.3p ). we 
explicitly show that there is a one-parameter set of solutions (I = 0, 1, 2, • • • , k — 1) in 
the "unitary" series q = p which can be tersely expressed as 

x = ^J^P ( rf'- 2J ^\ z ) yJ(z-l) l (z + l) k -\ (1.5) 

q f 2l-k 2l-k \ I T ] 7 

Q = pKezp^ ^ ' fc \z) y(z - l) k -\z + 1)', (1.6) 

in terms of the Jacobi polynomial^ P^^\z). 

This ansatz moreover enables us to explicitly describe the geometry of off-critical 
amplitudes. For instance, our ansatz for the (p,q;k) = (1,1; 3) case gives the following 
algebraic curve: 

F(x, Q) = Q 3 - 3/ixQ - x 3 = 0, (1.7) 

the k(> l)-cut cases, there are two kinds of indexes: (p, q) and (J>, q). The former corresponds to 
the standard CFT labeling of critical points; The latter corresponds to order of differential operators P 
and Q. These labellings are different in general [33] [19]. Since the Liouville continuum formulation of 
Z/c symmetric critical points has not been known so far, we only use (p, q) labeling except for the bosonic 
(k = 1) cases. 

2 There are two solutions in this system, that is, different phases (above and below the tip of a critical 
point) in the matrix models [17]. One solution is essentially the same as the bosonic case (|1.2[l with 
(p, q) = (p, q), which is called one-cut solution. The other is essentially different from the one-cut solution 
which is called two-cut solution (|1.3p . 

3 Here we neglect normalization factors and phases in front of polynomials which are irrelevant in this 
section. 

4 The definition and basic properties of the Jacobi polynomials are listed in Appendix [Dl 
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which contains precisely three symmetric cuts corresponding to the eigenvalue distribu- 
tions in the three-cut two matrix model. 

The organization of this paper is following: Some basics of the multi-cut two-matrix 
models are reviewed in section 12. 1L The scaling ansatz for the orthogonal polynomials 
and its relation to the u; 1//2 -rotated/real potentials are discussed in section 12.31 The 
solutions for critical resolvents and potentials are studied in section [2~4l (some examples of 
critical potentials are listed in Appendix [B]). The Daul-Kazakov-Kostov prescription [27] 
is reviewed in section 13.11 and the two-cut cases are discussed in section 13.21 The Z*. 
symmetric fc-cut cases are studied in section 13.31 Its algebraic structure is also discussed 
in section 13.51 Section @] is devoted to conclusion and discussion. 

Some of the results for the critical potentials are listed in Appendix [B] Direct evalu- 
ation of Douglas equation in the three-cut cases is given in Appendix O Several useful 
formulae about Jacobi polynomials are collected in Appendix [Ql 



2 Critical potentials in multi-cut critical points 

2.1 Preliminary of the multi-cut two-matrix models 

We start with a brief review of some elementary facts about the orthogonal polynomial 
system in the multi-cut matrix models, which also serves to fix the notation used in this 
paper. 

As is the usual two- matrix model [46], the multi-cut two-matrix models are charac- 
terized by the orthonormal polynomial system, 

a n {x) = -L (x n + ■ ■ • ) , (3 n (y) = -L (y n + ■ ■ • ) , (2.1) 

with an inner product, 

(a n (x),p m (y))= [ dxdye- Nw ^ a n (x) f3 rn (y) = 8 n>m , (2.2) 



defined with respect to the two-matrix potential w(x,y) = V±(x) + V 2 (y) — 9xy. We will 
use the normalization of the potential, 9 = 1, throughout the main text. The degrees of 
the potentials are chosen to be multiples of k, 

deg Vi(x) = kmi, deg V 2 (y) = km 2 . (2.3) 

Here k is the number of cuts and mi,m 2 are arbitrary positive integers. A distinct 
feature of the multi-cut matrix models is that the integration domain C x x C y is defined 
as a contour along radial directions specified by the Zk symmetry of the system in the 
complex plane [18], 

fc-i 

C x {=C y ) = |J w n Rc C, c/ = l, (2.4) 

n=0 

where u is a root of unity, which is chosen to be 

u = e— . (2.5) 
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The fc-cut matrix models possess a Z& charge conjugation, which is a generalization 
of the reflection transformation for the matrix (M — > —M) in the two-cut phases of the 
one-matrix models [15-17,37]. The Zj. charge conjugation in the k-cut two-matrix models 
is given [33] as 

(x,y) - {Jx^y) (j = 0,1,2,--- ,fc-l), (2.6) 

which preserves the interaction term, xy, in the two-matrix potential. In general, if the 
potential w(x,y) is a Z^ symmetric potential, 

w(u J x,uj~ 3 y) = w(x,y) (j e Z), (2.7) 

the system is said to be Z^ symmetric. The Z^ charge conjugation naturally defines a 
charge conjugation of the ortho normal polynomial as 

a^(x) = uT^>V), /#>(y) ee W n J'a n (w- J 'y), (j = 0, 1, • • • , k - 1). (2.8) 

In the Zfc symmetric cases, the orthonormal polynomials are invariant under this charge 
conjugation, equivalently we have 

a n (u j x) = u nj a n (x), (3 n (uj~ J y) = uJ- nj (3 n (y), (j e Z). (2.9) 

That is, the orthonormal polynomials are eigenf unctions of this Z& transformation. 

The orthonormal polynomial system is characterized by the recursive relation, which 
can be written as 

d 

xa n (x) = A(n,e- 9 ") ■ a n (x), N" 1 — a n (x) = B(n, e~ dn ) ■ a n (x), (2.10) 

ox 

y(3 n (y) = C(n,e- 9 -)-f3 n (y), N' 1 ^- f3 n (y) = D{n, e~ 9 «) • /3 n (y), (2.11) 

dy 



where 



fem 2 — 1 1 (Ami — l)(fcm2— 1) 



A(n, Z) ee y/R(n+l) V A s (n) B(n, Z) = -±= V £ s (n) Z s 

v -R(^ 



8=-l V "4 n J s= l 



(2.12) 

fcmi— 1 1 (fcmi— l)(fema— 1) 

C(n, Z) ee V^Tl) V C a (n) Z s , D(n, Z) = —== V D s (n) 

(2.13) 



8=-l 



and R(n) = h n jh n -\. Here Z = e dn is the index-shift operator [6] [47] [27] that shifts 
every index n by one unit on its right-hand side: 

Z m ■ a n+i (x) = a n _ m+i (x), Z l oA s {n)=A s {n-l)oZ\ (2.14) 

and Z T = Z~ l . Note that, by construction, A_i(n) = C_i(n) = (N/n) B\(n) = 
(N/n) Di(n) = 1 and R(n) cannot be zero, 

R(n)^0 (n = 0,1,2,- •• ,oo), (2.15) 
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including n — > oo. In the Z& symmetric case, invariance of the system asserts that the 
expansion coefficients A s (n), B s (n), C s (n) and D s (n) must be zero except for 

A mk _i(n), B mk+ i(n), C mfe _i(n), D mk+1 (n), (2.16) 

with m G Z and n G Z. 

From the recursive relation (I2.10p and (I2.1ip . one can check that the operator pairs 
(A, B) and (C, D) satisfy the canonical commutation relations, 

[A,B} = [C,D}= N- 1 . (2.17) 

Furthermore, taking inner product (12.21) between d x a n (x) and (3 m {y) and then integration 
by part leads to the following identity, 

B(n,Z) = V{(A{n,Z)) -C T (n,Z), D(n, Z) = V^(C(n, Z)) - A J \n, Z). (2.18) 

Combine these equations together, it can be seen that there is only one independent 
canonical commutation relation, 

[A,(-Cf r )] = [C,(-A T )]=N- 1 . (2.19) 

One can also extract equations for A and C from eq. (12. 18ft : 

A(n, Z) = k(C T (n, Z))} - Z~'-^= (2.20) 
C(n, Z) = \v{(A T (n, Z))} - Z^^= (2.21) 
Here the bracket is defined as 

The system of equations was originally used to extract the string equations [4, 5] and 
can be considered as the fundamental equations dictating the double scaling limit of the 
two-matrix models. Once we solve A and C, eq. (I2.18P leads to the solution of B and D 
in straightforward manner. 

The strategy for solving critical points in two-matrix models was proposed in [30,31]. 
The observation is this; the operators A and C have the following scaling behaviors at 
(p, q) critical points (determined by critical potentials): 

A(n, Z) ~ a p/2 d p + • • • , C T (n,Z) ~ a q/2 d q + ■ ■ ■ , (a -> 0), (2.23) 

where the lattice spacing a is introduced to satisfy 

N =g str a2 ->0, — = exp(-ta 2 ) -> 1, 

9 n = -a 1/2 g str d t = -a 1/2 d -> 0, (a -> 0), (2.24) 

and the operator 9 n is the index-derivative operator in the index-shift operator Z = e _s ™. 
Conversely, if one assumes this critical behavior (I2.23p . then there exist critical points. 
That is, we can obtain critical potentials which satisfy eq's (12.201) and (12.211) [30,31]. 
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Explicit solutions to this problem have been given in [27] for the general (p, q) critical 
points in the one-cut cases as follows: For each pair of potentials with degrees (mi, m 2 ) 
(k — 1 of eq. (12.31) ). the maximal order of the differential operators (12.231) is (p, q) = 
(m2,mi) and the operator A and C is given as 

— (1 - Z) m2 , — (1-Z) mi 
A(n, Z) - A(*, Z) = ^/Rj C(n, Z) - C(*, Z) = ^R^- ^— , (2.25) 

in the large N limit. Here R* is the critical value of R(n), 

R{n) -> i?*, (n -> oo). (2.26) 

Then it can be shown that, with these critical forms of A and C, we can solve for the 
critical potential w(x,y) which satisfies eq. (12.201) and (12.211) . 

In the rest of this section, we generalize this idea to the multi-cut setting. 



2.2 Smooth functions in multi-cut critical points 



The procedure of section 12.11 implies that the orthonormal polynomials become smooth 
functions of the index n at the vicinity of the critical points, especially at eq. (I2.24p . 

The scaling ansatz for the multi-cut cases is to find a sequence of orthonormal polyno- 
mials ot n (x) (and /3 n (y)) which are smooth continuous functions in not only x but also the 
index n. First of all, since the orthonormal polynomials behave under Z& transformation 
as 



a n (x 



a n {ujx) = L) n a* n (x) ~ u n a n (x) + 



(2.27) 



there should be k distinct smooth scaling functions, each of which has a distinct eigenvalue 
of the Zfc transformation. This is obvious in the symmetric background (12.91) . Hence, 
we separate the polynomials into the following k sets of polynomials 



{a kl+j (x)} leZ (j = 0, 1, • • • , k - 1), 



(2.28) 



with respect to mod k. For each set labeled by j, we assume that there exists a smooth 
function in I which can be regarded as the continuum limit of the sequence. The proper 
choice of the smooth functions was first proposed in [18] and we further write it as 



i>i(x) = 



( ipi,o( x ) \ 

y^ jfc _i(x) J 



xi(y) = 



/ xi,o(y) \ 
xi,i(y) 

\xi,k-i(y)J 



(/ e Z), 



(2.29) 



with 



kl-\-j hi | j 

uj —a M+j (x), xi,j(y) =uj~ p k i+j(y)- 



(2.30) 



The phase factors indicate that the polynomials are approximated by smooth functions 
in the directions of x = uj m+1 / 2 x (— > oo m+1 ^ 2 oo) and y = uj~( m+1 l 2 ^ y (— > cu~^ m+1//2 )oo): 



m+l/2~ 



LO 



X) ~ 



:X kl+j + - 



Xi,j( 



LO 



-(m+l/2)~ 



LO 



-mj 



ykl+J + 



(2.31) 
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since the phase does not depend on / here. These directions are actually the Liouville 
direction [28,29] [37] [33]. The main reason for this ansatz is that the critical point of 
eigenvalues, (x*,y*), is at the origin = (0,0) when we consider Z fc symmetric 

potentials w{uox,uo~ l y) = w(x,y)U 

The consequence is that the operators acting on the vector- valued polynomials ipi(x) 
(and Xi(y)) become smooth functions of /. The operators in this basis are expressed as 



xi/> l (x)=uV 2 A(l,Z)-il> l (x), 



N- 1 ^ l {x)=u-^ 2 B(l 1 Z)-^ l (x), 



y X i{y) = ^ 1/2 C{1, Z) ■ Xl (y), iV" 1 ^ Xl (y) = u l / 2 V(l, Z) ■ Xl (y), 
with explicit expressions: 



jfc-i 



A(l,Z) 
B(l,Z) 



Km+iV] diag ( V A s (kl+j)uj-^ Z s ) xT(Z) 
r=0 o<j<k-i\ ekZ+r t 

1 fe-i 

= J2 diag ( B a (kl+j)u-^Z') xT(Z) k - 
fe-i 

c(/,z) = v / ^^E dia £ ( E c a (ki + j)u^z')xr(z) k - 

r=0 — J — s£fcZ+r 
fc-1 

V(l,Z) = ^==J2 diag (J2 D s {kl+])u^ Z°)xT(Z) k ~ r 
V rCkl r, 0<i<k-l \_Tr, _ 7 



fc-r 



where 



with 



'n 



kl+l 



diag (v^RpTTTT; 



diag 



VTZm 0<i<fc- A a/R(A^ + j 



Mo 



diag ( A 

o<i<fe-i 



Ai 



A 



it-i / 



and 



r(z) = 



/ o 



i 



z 



z 





r(z) fc = 4 



Here Ik is the identity matrix of k x fc. Now eq. fl 2 . 1 8 [) becomes 

B = u 1 / 2 V± L X ' 2 A\ - C T , 2? = V? (^" 1 / 2 C) - A T , 

One can explicitly check this fact from the critical potentials we list in Appendix iBl 



(2.32) 
(2.33) 

(2.34) 
(2.35) 
(2.36) 
(2.37) 

(2.38) 



(2.39) 



(2.40) 



(2.41) 
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in terms of the k x k matrix- valued operators A, £>, C and T>. 

Note that the smoothness of polynomials ipi(x) and Xi{v) i n the index I only implies 
smoothness of sequence of recursion coefficients {A s (kl+j)}il Q in the index I. This means 
that each recursion coefficient sequence (labeled by j) can be approximated by a smooth 
function with the following Taylor expansion with respect to the scaling parameter aa; 

A s (kl + j) = A s (k^+j) + a 1 / 2 f«\t) +ag<?\t) + ■■■. (2.42) 

Here the lattice spacing a and scaling parameter t are defined as 

N =g str a2 ->0, — = exp(-ta 2 ) _> l ; 

^ = -a^g^dt = -a l ' 2 d 0, (a 0), (2.43) 
for the (p, g) critical points of the multi-cut matrix models. 

2.3 The cj^-rotated v.s. real potentials 

It has been argued in the two-cut matrix models [14] that Z 2 symmetry breaking critical 
points can only be realized in the potentials V(x) which has an imaginary number i in 
the following way: 

V(x) = V (-ix). (2.44) 

Here the potential Vq(x) is polynomial in x with real coefficients. This consideration 
has a natural correspondence in the Liouville theory because the Liouville boundary 
cosmological constant ( is related to the eigenvalue x with an imaginary number i [17] as 

x = zC. (2.45) 

A natural extension of this consideration is the a/^-rotated potential (12.4=6)) . From our 
resultsj§ within the prescription reviewed in section 12.11 one can a ls° find real-potential 
solutions for the Z2 symmetry breaking critical points. Therefore, we consider both cases 
in this section. We will show in the next subsection that, at least in the sense of critical 
potentials, the map between two potentials is just an analytic continuation. 

2.3.1 The cj 1 / 2 -rotated-potential models 

The cj 1 / 2 -rotated potentials are defined as 

w(x,y) = w Q {u- 1 ' 2 x,uj 1 l 2 y), (2.46) 

with a real-coefficient potential Wo(x,y). If one puts k = 2, then this turns out to be 
eq. (12.441) given in [14]. From the definition, this potential satisfies 

w*(x,y) = w (lux, (2.47) 

6 See Appendix [H 



S 



and consequently one can show thaf] 

S m , n = (/3 m {y),a n (x)) = (/3 m (y) , a n (x)}* 

dx*dy* e -^W) p* m ( y *) a * n ( x *) = f dxdy e'^^ f3* m (y) a* n (x) 

CxC JCxC 

dxdye- Nw ^ P* m {uy)a* n {uj- l x) 

ICxC 

= (u- m {3* m (uy),u; n a:(u~ 1 x)), (2.48) 

therefore the uniqueness of the orthonormal polynomials implies 

a* n (x) = u- n a n (ux), (3* n {y) = u^u^y). (2.49) 

Here we have used the invariance of contour C* = C and the invariance of the mea- 
sure, d(ujx)d(uj~ 1 y) = dxdy. This complex conjugation corresponds to the Z& charge 
conjugation of the polynomials ( \2.8h . 

<(x)=a«(x), (3* n (y) = ^)(y). (2.50) 

The hermiticity of coefficients in the operators A, B, C and D is also obtained in the 
same way as 

A,(n)* = uj~ s ~ x A s {n) <^> Ajn) = uj~^~ A s {n) £ R, 
BJn)* = ou- s+1 BJn) BJn) = BJn) G R, 



r- x A a (n) 




A s (n 


r s+1 B s (n) 




B s (n 






C s (n] 


u'- x D t {n) 




D s {n 



D s (n)* = uj s ~ l D s (n) D s (n) = uj^ D s (n) E R. (2.51) 

Comparing this with the expression (12.341) , (I2.35P , (12.361) and (I2.37P , one can see that the 
coefficients of the matrix-valued operators A, B, C, T> are all real number. In this sense, 
eq. (I2.4ip can have critical solutions with this d// 2 -rotated potential, including the Z*. 
symmetry breaking critical points. 

2.3.2 The real-potential models 

The real-potential models are defined by real two- matrix potentials w(x,y) with all co- 
efficients real: 

w(x,y)eR, (x,yeR). (2.52) 

The same discussion as above shows that the polynomials a n (x) and /3 n (y), and the 
operators A, B, C, and D are all real functions. 

The relation between eq. (I2.18P and the smooth operators becomes simpler if one 
considers another convention of smooth orthonormal polynomials]^ 

^r\x) = (-i)^cw*), xtP l \y) = (-i^/Wv), (2-53) 



It can be similarly shown that h n are real functions. 

8 Thc two-cut matrix models are studied almost with this convention in literature. 
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instead of eq. (12.301) . The corresponding smooth operators A^ Tcal \ J g( real ) ) C( real ) and 2)( real ) 
are then given as 



(2.54) 



yx { r\y) = c (rcal) (^, z) ■ x [ r\y), n^-^- x f cal) ( y ) = v^\i, z) • x { r\y)- 

(2.55) 

Here ipf cal \x) and Xi (v) are the counterpart of ^(x) and given in eq. (I2.29p . 

One can easily see that the operator v 4( real ) can be obtained from A by just replacing the 
phase t ) ;~( s + 1 )/ 2 in eq. (I2.34p with (— l)H s+1 )/ fe ] (the phases in the other operators B,C,V 
are also similarly replaced: the phase cu / 2 is replaced by (— l)[ / fc ]) and by replacing the 
matrix T(Z) with the matrix r( real )(Z), 



/ Z 

Z 



T{Z) 



\ 



\ Z k-i 



z 





r (rcal) (z) 



/ z 

z 



-1 
\z k ~ l 



\ 



z 





/ 

(2.56) 



which satisfies (r^ real ^(Z)) fc = — 1^. Here [— (s+l)/k] is the Gauss symbol, which expresses 
the maximal integer less than or equal to — (s + l)/k. Then eq. (12.411) is equivalent to 



g (real) = ^/^(real)) _ Q 



(real) - 



real) - 



(2.57) 



Since everything is real here, one can have critical points (even symmetry breaking 
ones) in the same sense as the a/^-rotated models. 



2.4 Multi-cut critical potentials and resolvents 



The large N solutions to eq. (12.411) are given by critical values of the matrix-valued 
operators A and C (defined in eq. (I2.34p and eq. (I2.36P ). 



A(n, Z) A(*, Z) ~ dP 



C T (n,Z) ^C T (*,Z) ~ 6P 



(2.58) 



and the critical potentials w*(x,y). At the vicinity of critical points, the coefficients of 
operators A (and C) should be approximated by the following smooth functions; 



A s (kl+j) = A s (koc+j) + a 1 ' 2 fU\t) +ag?>(t) + 



(I 



oo 



(2.59) 



In principle, the leading coefficients A s (koo+j) (and C s (koo+j)) can depend on the index 
j and the critical operators of A (and C T ) are expanded by these coefficients. For a while, 
let us see what happens in some special cases@ That is, we consider the (p, q) = (m 2 , mi) 
critical points realized in the potentials of degree (kmi, km^) (See eq. (12. 3p ). 



3 The calculations in this section have been carried out mostly with the help of Mathematica™ . 
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First we solve the condition (12.581) . Since A and C T are almost similar, the operator 
A is only considered here. By solving eq. (12.581) in this critical point (p, q) = (m 2 , mi)F°l 
one can obtain the critical operator of A as 

A(*,Z) = y/n^ ^ l ~z )P T{Z). (2.60) 

The matrix y/Hkoo+i is defined in eq. ( 12.381) and the elements R(koo + 1 + j) can depend 
on the index j in principle. From eq. (I2.34p . therefore, one can obtain 

I £ Mkoo+j)uj-^Z^= { ^-fl (2.61) 
sekZ-i 

Here each term contains a phase factor . Since s+1 is a multiple of k (say s+1 = kl) 
and a; - ^ = e ml = (—1)', we get 

y £ A s (koo+j)Y s+1 = {1+ ^ k)P (Y k = -Z k ). (2.62) 

s+iekz 

That is, the operator A (and similarly the operator C) is given as 

A(koo + j, Z) = y/R(koo + l+j) ( 1 + Z k )\ 

C(koo + j, Z) = y/R(koo + 1 + j) (1 + Z " )9 . (2.63) 

Z 

We then consider eq. (12.201) and eq. ( 12.211) . The above operators A(koo+j) and C(koo+j) 
(j — 0, 1, • ■ ■ , k— 1) satisfy eq. (I2.20p and eq. (12.211) for each j. By solving these equations, 
one can obtain the critical potentials and the solutions are realized when one has 

R(koo + 1) = R(koo + 2) = • ■ ■ = R(koo + k) = (2.64) 

Therefore, the critical behavior of matrix-valued operators A and C is finally given as 



A{*, Z) = y/R S z ' T(Z), C v (*, Z) = ' r(Z) k ~\ (2.65) 

Several examples of critical potentials calculated by using the above procedure are listed 
in Appendix [B] and some critical potentials are drawn in Fig. [TJ Note that all the 
parameters here, including i?*, are fixed in this procedure0 



10 By using the expression Z = e dn , one performs the Taylor expansion of A(*, Z) in terms of d n 
and requires A(*,Z) = (const.) x df L + 0(df l +1 ). Then one can obtain constraints on the coefficients 
A s (koo The result is eq. (|21)0|l 

n The counting of redundant degree of freedom is the same as in the one-cut cases [27]. It can be 
seen as follows: Originally we have four redundant degrees of freedom, (x, y) — > (ax + b,cy + d), but the 
critical points of eigenvalues are located at the origin, (x*,y*) = (0,0), we have chosen the scale, 0=1, 
and (x,y) — * (ax,a~ 1 y) was fixed by choosing the convention of overall coefficient in the orthonormal 
polynomials (|2.1[) . In this sense, there is no free parameter here. One can also say that one can freely 
tune the parameter by turning on 9. This procedure is a redundant deformation which does not 
change the result, and we will use this fact in Appendix [Bl 
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The result (12.641) implies that all the k coefficient sequences A s (koo+j) (j = 0, 1, • • • , k- 
1) have the same asymptotic behavior at a — > limit, i.e. at critical point, the coefficients 
satisfy 

A,(koo + j)=A B (*), C s {koo+j) =C 8 (*) (2.66) 

even though A s (n) can be approximated by k different scaling functions (e.g. f^\t)) 
in the sub- leading terms of eq. (12.591) . This was just an assumption in the previous 
discussion of the two-cut two-matrix models [33] but here we have checked that this is 
generally correct in the multi-cut two-matrix models. 

Next consider the critical points of (p,q) = (m2,mi — 1) realized in the potentials of 
degree (kmi, km^) (See (12.31) ). This contains symmetry breaking critical points (p, q). 
The critical operators of A and C are given as 

A(t, Z) = - z r r(Z), (2.67) 



cT(t Z) = ^(i-nii-g^W? , r(z) .- 1+ 

Zj 



%J2°r(*) {1 ^ nzy, (2.68) 



r=0 



with C r (*) = uj 2 C r (*) (see eq. (12.511) ) . In particular, from eq. (12.36j) . one can see the 
following relations: 



Z ' Z 



s+lGfcZ 

u r -^Z r C s {*)uo^r Z s ~ r = C r {*){\- Z k fZ r . (2.69) 

s—rdkli 

Therefore, the operators A and C are also obtained in the similar manner as eq. (I2.63p . 
and given as 

C(», Z) = ^ ft + ^'U + ^-iW^) + g Cr (») . (l + Z')V. (2.70) 

r=0 

The critical points and potentials in this cases are also obtained with these operators 
and eq's (12.201) and (12.211) . which are also listed in Appendix [Bj Note that the critical 
potentials now depend on the parameters {CV(*)}^~q and Ck m2 -i{*)- 

Note that the difference between the u; 1//2 -rotated potentials (12.461) and the real po- 
tentials (12.52j) is almost the complex phase of C r (*) = C r (*) e m(r ) (CV(*) g K; r = 
0, 1, • • ■ , k — 2) in eq. (12.701) . Since we obtain the critical potentials by solving eq's (12.201) 
and (12.211) by using eq. (12.701) . the map from real critical potentials to k/^-rotated ones is 
also given by just analytically continuing the parameters {C r {*)}^Z^ in the real potentials 
to the complex values {C r (*) ^~ ( - r+1 ^ 2 }r=o- This means that all the symmetric (p,q) 
critical points and Z^-symmetry breaking (p, q) critical points can be realized within 
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(a1) 



(a2) 





(b2) 



(b2) 





(d) 



(c2) 



Figure 1: The critical potentials of (a) (p, q; k) = (2, 2; 3), Vi(a 



V 2 (x) 



~3~ 



6 

10x 6 



; (b) (p,q;k) 



(al), 



(2, 2; 4), Vi(ar) = V 2 (x) = -3a; 4 + — ; (c) (p,q;k) = (2, 2; 6), V^x) = V 2 (x) 

(bl) and (cl) are plots in the real axes. The figures (a2), (b2) and (c2) are plots of the real part of the 
potentials Rc(Vi(x)) on the complex sheet of x. The straight lines with arrows on them correspond to 
the contour C x of integration. Darker color means lower height. 
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both real potentials and cu^-rotated potentials of fc-cut two-matrix models, and also 
that these real/a; 1/ ' 2 -rotated critical potentials are related to each other by just the above 
analytic continuation. 

One can then realize the canonical operators in terms of the Lax pair, and they are 
given by 



A(n, Z) = a^ /2 P(t, d), C T (n, Z) = -a« /2 Q(t, d), 



XUJ 



-1/2 _ _p/2 



with 



p-i 



p(t, d) = v^Ki-ky x r + H n(t) 



a p/2 C (2.71) 



(2.72) 



n=0 



k-2 



Q(t,d) 



(i + M*))r w + JaWr r dt + Y,H n (t)d«- n , (2.73) 



9-1 



r=0 



n=0 



and the canonical commutation relation (1 2 . 1 9 [) becomes the Douglas equation, [P, Q] = 
g s t r Ik- The functions H n (t) and H n (t) are some combination of scaling functions in f)2.59p . 
The parameter a is a lattice spacing of d n = — a 1 / 2 ^,.^. Here we define 



r = nr 



/o 1 

1 



V 



\ 



1 



r fc = h. 



(2.74) 



Here H n (t) and H n {t) G 1R are all real functions. Since we can erase some parts of 
sub- leading scaling functions H n {t) and H n (t) by changing normalization of ortho normal 
polynomials (some similarity transformation), we can always put 



'5 -M 5 



i=l 

k — 1. Here we extend the meaning of indices as [i?i(t)] . +lk 



j+mk 



for / = 0, 1, 

It is also convenient to use the diagonalized basis of T which is useful for analyzing 
multi-component KP hierarchy [33]: 



/I 



Q = U ] TU 



\ 
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\ 



(2.76) 



uj k - 1 J 



and the matrix elements of U are given by 



Ui 



10 



"71 ' 



(2.77) 
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In this basis, the relation (12.751) is simply given as 



[rfH!{t)U] =0 (i = 1, 2, - - - ,k). (2.78) 



However, from the analysis we have shown in section I2.3[ the component functions in the 
operators ^( KP ), £( KP ), C( KP ), £>( KP ) ^(kp) = tfxU) are complex-valued functions now. 

Here we comment on the "unitary" series which is defined by the condition of V\(%) = 
Vr>(x)r^ We should note that the unitary series of the symmetric background is given 
by (p, q) = ip,p) instead of (p,p+l). This can be understood as follows: In general, the 
"unitary" condition gives 

A(*, Z) = v 7 ^ 1 z ] r(Z), C T (*, Z) = z _ x > Y{Zf-\ (2.79) 

in the Z& symmetric background. Here m = mi = m 2 . In the superstring {k = 2) cases, 
since the matrix T(Z) satisfies 

Y{Z) k ~ l = T(Z), (2.80) 
the following combination of the operators A and C T 

i (A(*, Z) + (-l) m C T (*, Z)) = a m x y/R~*(-k) m r «9 m + • • • , 

^ (U(*, Z) - (-l) m C T (*, Z)) = -a m+1 x y^(-k) m+1 T «9 m+1 + • • ■ , (2.81) 

gives the (p, q) = (m,m + 1) system. The same argument also applies to the bosonic 
(k = 1) cases. For the k (> 2)-cut Z fc symmetric background case, the matrices T(Z) and 
T(Z) -1 are essentially different matrices; hence, the previous argument fails. This means 
that the unitary series should be (p, q) = (m, m). 

On the other hand, if we consider the Z& symmetry breaking critical points which 
correspond to minimal fractional superstring theory [19], 

A(*, Z) = y/R ^ l ~^ Y{Z) =a?x ^/W*(-kf T & + ■ ■ ■ , 

Zj 

c T (*,z) = v^ (1 | ] T(Z) + v^ V Z J nz) k ~ l 

= afi x y/R^tfrtf + ■■■ , (2.82) 

the matrices of the leading terms in A and C T are both given by T(Z). That is, the 
situation is similar to the super(bosonic)string cases. In fact, unitary minimal fractional 
superstring theory corresponds to the (p,q) = (m,m + 1) critical points, instead of 
ipA) = (m,m). 



12 This is also called "Z2 symmetric" case in literature, in the sense of exchanging the matrices X «-> Y . 
Since the terminology of "Z2 symmetric" is dedicated to symmetry in this paper, we will use the term 
"unitary" referring to the V\ (x) = V%(x) case. Of course, the terminology of "unitary" should come from 
the worldshect CFT of dual string theory, and therefore our terminology "unitary" should be understood 
as a formal notion. In the case of one- and two-cut matrix models, this situation always gives unitary 
sequences although obtaining unitary sequences does not necessarily require this condition. 
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3 Macroscopic loop amplitudes 



In this section, we consider off-critical behavior of the critical points we identified in the 
multi-cut matrix models. The prescription we adopt for studying off-critical amplitudes 
was developed in the one-cut two-matrix models [27], which we refer to as the DKK 
prescription. In the following, we first review their procedure then consider the extension 
to the multi-cut cases. 



3.1 The one-cut cases: a review of [27] 

At the (p, q) one-cut critical points, the canonical pair of A and C operators was given 

as 

Here we choose = 1. Considering the dependence on n/N (= exp(— ta 2 ^ )) in the 
operators A and C is equivalent to turning on the coupling of the most relevant operator, 
t. Since we have dimensionful parameter t, we can make a dimensionless combination 
with dt as 

z = Ktdt, [z,t n }=nnt n , (3.2) 

p+q 

with k = g s tr/t p+q - 1 ■ Then the off-critical amplitudes of the canonical operators are 
generally expressed as 

A(n,Z) = a* \ p tt p (z), -C T (n, Z) = J\ q £ q (z), X = t*£=*. (3.3) 
Note that the leading behavior of the polynomials n p and £ q are 

7r p (z) = (-If z*> + ■ ■ ■ , i q {z) = -z« + ■ ■ ■ . (3.4) 
So the commutation relation becomes 

N- 1 (= K a-¥\ p+q ) = [C T ,A] =N- i q A ^ _ * , (3.5) 

and 

qn' p (z)^ q (z) -pZ q (z)n p (z) = (p + q-1). (3.6) 

This was found in [27, 42] and is one of the main equation for off-critical amplitudes in 
the two-matrix models. 

It is known that this system has a simple solution for conformal background which is 
given by the Chebyshev polynomials of the first kind, T p (cosh(r)) = cosh(pr) [22]. In the 
case of unitary series (p, q) = (p,p+ 1), the equation (13.61) becomes the addition formula 
for the hyperbolic sine and cosine functions: 

m , , Nm , . mi / \ r-r-r i \ shah ( vt ) cosh ( or ) — sinh ( vt ) cosh ( or ) 
qT' p {z)T q {z) - P T' q {z)T p {z) = pq ^ Lhr 

sinh(g — p)r 

= -P<1 ^-r^ = -pq {q=p + l). (3.7) 
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Consequently, the solution for q = p + 1 is 

A(n, Z) = at x 2 (-^X^ T p {z / c) = a% x 2 (-^) P cosh(pr), (3.8) 
C T (n,Z) = a^ x2(£\) 9 T,(*/c) = x 2 (|a) 9 cosh(gr), (3.9) 

with the normalization factor c which satisfies 2(—l) p (c/2) p+q ~ 1 = (p + q — l)/pq and 
^ = ccosh(r). 

For the general (p, q) cases, the Chebyshev-polynomial is not a solution to eq. fl3.6j> . 
However, aside from the unitary series, the most relevant operator and cosmological 
constant are different in general; considering this fact here, it means that there can be 
a Jacobian factor in the eq. (13.61) in general. Therefore, we will utilize this degree of 
freedom to deform the eq. (13.61) ; it was proposed in [27] that the simple solutions with 
the Chebyshev polynomials can be derived from 



rr> ( \rr f \ T I \ r r l i \ sin(gr) cos(pr) - sin(pr) cos(gr) 

V T i {w)T^L- (w) - q T_jl (w)T p (w) = pq —— r 

q -p q-p q-p q-p — pj sm{q — p)T 

pq 



q-p 



(3.10) 



Here T_s_{w) = T q (z(w)) with w = cos(g — p)r = T q _ p (z). This procedure corresponds 
to utilizing the canonical pair (w,/i) of the cosmological constant fi [27], 

[W7,/I n ] = KUf/ 1 , (3.11) 

instead of the canonical coordinates (z, t) of the most relevant operator t, 

[z,t n ] = nnt n . (3.12) 
So in general the equation we have is 

i n' q {z) C p (z) - 1 n q (z) t' p (z) = w' p>q (z). (3.13) 

In this procedure, the requirement of matrix model is that the functions of z in the 
operators A and C T are all polynomials in z. One can expect that w Ptq (z) should be 
polynomials in general background. 

3.2 The two-cut cases: revisited 

The solution for conformal background in the superstring cases was first obtained in [34]. 
There are two classes of the solutions. The first class is called the one-cut-phase solution: 

C = v^cosh(pr) = y/jiTpiz), d ( Z(C(z)) = ^ cosh(gr) = ^ T q (z). (3.14) 
The second is called the two-cut-phase solution: 

C — v^/isinhQoT) = y/jlU p (z)V z 2 — 1, d^Z(((z)) = fi 2 ? sinh(gr) = fi 2 p U q (z)v z 2 — 1, 

(3.15) 
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Here we use z = cosh(r) and Z(Q is the FZZT disk amplitude. Realization of this 
solution in the pure-supergravity critical point was shown in [17]. Realization of (p,q) 
solutions in the two-cut two-matrix models was shown in [33], with the SFT formulation 
[48] [28,43] and its prescription utilizing the Wi +00 constraints [32]. 

In the DKK prescription, what we need to consider is the Douglas equation for 
the matrix- valued Lax operators A(n, Z) and C T (n, Z) of eq's (12.341) and (12.361) . Since 
A(n, Z) and C T (n, Z) can generally be expanded in terms of Pauli matrices <7j, one can 
decompose the string equation in the following form: 

= [C T ,A] a-^X-W* - nl 2 = (CCRO) k° + (CCR1) «J 2 + (CCR2) na 3 + 0(k 2 ). 

(3.16) 

The one-cut solution is trivially the solution found in the bosonic case. The two-cut 
solution is less trivial but can be realized as 



A=aV 2 \c- 1 (-cXy,U^ 1 (z/c) 



z + c 
z-c 



C T = a«/ a {c- 1 ( C A) tf ,l7 4 _ 1 (z/ C )^° c ^ C )}, (3.17) 

with 4(— l)^c^ + ^ _1 = (p + q — l)/pq. One can easily show that this satisfies the Douglas 
equations (I3.16p . The curly bracket {, } denotes anti-commutator which is necessary to 
satisfy (CCR2). 

The implications of this system are the following: 

• The zero-th order equation of the Douglas equation, (CCRO), indicates that A 
commutes with C T . That is, A and C T can be diagonalized simultaneously!]!] 

• Eigenvalues of the matrices A and C T satisfy the same equation (13.61) as in the 
bosonic case. The important difference is now the eigenvalues can include non- 
polynomial parts (e.g. the square root in eq. (13.151) ). as long as all the matrix 
elements remain to be polynomials in the matrix-model realization (13.17p . 

• Usually, some naive expression of the matrix-model realization (e.g. eq. (13.171) 
without the anti-commutator) cannot satisfy (CCR2). In this case, we need to add 
proper first order corrections to the operators A and C T 



3.3 The /c-cut cases: symmetric background 

Now we can extend the previous discussion to the fc-cut matrix models. Here we will 
focus on the general considerations; more detailed investigation of the three-cut string 
equation can be found in Appendix O For the sake of simplicity, we consider the models 

13 It has been shown in [33] that each pair of eigenvalues, (Pi, Qi) (i = 1, 2), corresponds to the FZZT 
brane amplitude with each different R-R charge. 

14 In this case, all we need to do is just to change the ordering of A and z. But in the higher-cut cases, 
we need to add purely first order corrections to the scaling functions in the operators A and C T . See 
section ICl.lt 
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with the Zfc symmetric background. In this case, the operators A(n, Z) and C T (n, Z) are 
restricted to 



A{n, Z) =a i>/2 \ p 



( F x {z) 

F 2 (z) 



\F k {z) 

I o 

G 1 (z) 

G 2 (z) 



\ 



F k -i(z) 




(3.18) 



G k {z)\ 



\ 





G k ^(z) ) 



(3.19) 



The matrix-model requires that all the functions F^{£) and Gi(z) are polynomials in z. 
Thus the diagonalization is just given by 



A{n, Z) ~ a?' 2 \v n p (z) Q, C T (n, Z) ~ -a« /2 A« 



(3.20) 



with 



7T 



l/k 



l/k 



l/k 



and hp(z) = Yli^i( z ) anc ^ hq( z ) = TliGi(z). Here p p (z) , rjq(z) , u p (z) and are all 

polynomials and the equality "~" indicates diagonalization. Then eq. f)3.6p is expressed 
as 



q tip(z) hg(z) - p h' q (z) hp(z) =k(p + q-l) (hp(z) h 9 (z)) 



k-l 
k 



(3.22) 



Since the left-hand side of the equation is a polynomial, we need to resolve the fc-th root 
(...)( k ~ 1 ^ k in the right-hand side. A natural ansatz is 

up(z) = {z- bf-\z - c)\ v q (z) = {z- b)\z - c) k -\ (3.23) 

which is characterized by the two indices (k, I). Then the equation turns out to be 



(z -b)(z- c) qp'Jz) rj q (z) -pv'q(z) Pp(z) 



+ 



+ 



(q - p)z + pr/p - q-f q pp{z) r} q (z) =p + q-l, 



with 



7j5 



(k-l)b + lc 



Iq 



lb+(k-l)c 



k k 
This solution can be realized in the matrix model as 



A(n, Z) = aF /2 \Ppp(z) 



Yz-M 



C T (n,Z) = -a^ 2 \%{z) 



Y~ x z - M 



(3.24) 



(3.25) 



(3.26) 
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with 



/O c 




M = 



\ 



c 




b 

0/ 



/o 

b 



M = 



b 

c 



V 





c 0/ 



(3.27) 



Note that there is a sequence of / c's in M and a sequence of I 6's in M. These matrices 
satisfy M k ~ l = ^"'"V^M and [T,M] = [Y-\M] which result in (CCRO). Of course, 
there is some ambiguity of taking permutation among the positions of b and c0 Also we 
have the constraint of eq. (12.751) . which means that 



+ 



(3.28) 



With the above setting, we can solve equation (I3.24p . From straightforward evaluation 
of eq. (I3.24p . one observes that the parameters b and c satisfy 

b + c = 0, (3.29) 

except for p = 1 case; it can also be seen that 7^ is related to rjq(z) as 

% (z) = -A;«(^- 1 + 7 ^- 2 + .. 
For unitary series p — q, one can see that they satisfy 

Vp(z) = Pp(-z), 

which leads to simple solutions of eq. (I3.24p ; for instance, we list first five solutions of 
(k, I) = (3, 2) in the following^ 

p 2 W = (-3) 2 (, + ^), 



(3.30) 



(3.31) 



PaU) 

p4 z ) 



(-3) 3 U 2 + ^-^ 



(-3) 4 (V + ^ 



-3) 5 (z* 



„3 



")■ 

27/ 

5c^ 

9 Z 81 
17c 2 2 97c 3 



7c s 
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-z — 



567 



z + 



128c 4 
1701 



(3.32) 
(3.33) 
(3.34) 
(3.35) 



15 See Appendix [Cl where we explicitly study all possible solutions of the string equations in the three- 
cut cases. 

16 It should be stressed that the ansatz (|3.23|) can generate polynomials (i.e. solutions of eq. (|3.24p ) in 
all the (p, q) critical points not only in the unitary cases p = q listed below. 
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Our claim is that the solutions for the unitary cases (p = q) can be expressed by the 
Jacobi polynomials, pff* (z). In the diagonalized form (13.201) . the solution is given as 



A{n ' Z) ~ ^ 2cT(2p- 1) Wh^^) y(z-c) l {z + c)*-' x ft 



C T (n,Z) ~ a tf/ 2 (2cfeA) p Pj:r^'^^z/c) SJ(z-c) K -\z + C y x n- 1 , 

(3.36) 

with 

{-iy(2ck) p+q 1 = " ^ _ 1} — p^)^ + 2^, p - *¥)B(p, p) ■ (3 - 3?) 

Here B(x,y) = T(x)T(y)/T(x + j/) is the beta function and also ft is given in fl 2 . 76 [) . 
Note that the relation (I3.3ip is essentially the identity (ID. 71) of the Jacobi polynomials, 
Pn^\z) = (— l) n Pn' a \— z). We show that the above expressions of A and C solve 
eq. (13.241) for general indices (k, I) in Appendix ID. 41 

3.4 Some other solutions for the p = 1 cases 

For the three-cut critical points of (1, q), some of the solutions can be written in terms of 
the Chebyshev polynomials. Since this solution is very similar to the solutions of one-cut 
and two-cut cases, it provides simple examples to see what happens in multi-cut cases. 
In the case of p = 1, the constraints (I2.75P turns out to be 

7s= (^6 ± ic = (3 38) 

When we have (k, I) = (3, 2), then b + 2c = 0. A solution which satisfies this constraint 
is given by 

A(n,Z) ~ a 1/2 (-3)A(2cosh(^r)) i x ft = a 1 ' 2 (-3)\yf(z - c) 2 (z + 2c) x ft, 
C T (n,Z)~a^ 2 c- 1 (^)"(2cosh(^r)) 1 cosh(^-U) x ft" 1 

(it)* w ^ ( " /c) v(^- c )(^ + 2c ) 2 x q_1 > ( 3 - 39 ) 

with z = ccoshr. Here the polynomial W n (z) is the Chebyshev polynomial of the fourth 
kind: 

coshf 2n+1 r) 

H/ n (cosh(r)) = TTTV^' W M = T n{z) ~ {z - l)U n ^(z). (3.40) 

v ' cosh(^) 

This satisfies the deformed equation (I3.13p with the canonical pair (w, /j)f^l 

[w,fi n } = nnfi n , w{z) ~ cosh(g - 2)r ~ T\ 9 -. 2 \(z / c) , (3.41) 

with some proper normalization. 
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Note that there is no (1, 2) critical point in the Z 3 -symmetric 3-cut cases. 
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3.5 Geometry of the multi-cut solutions 

Here we consider geometry of the multi-cut solutions we have obtained. The purpose 
is two-folded; first, it serves as a check of our ansatz ( 13.231) ; second, it provides some 
physical intuition to the multi-cut system. The algebraic equation of the p — 1 solution 
(I3.39P is easily written down as 

F((, Q) = 2^C 3/2 (T 3 (Q/2^C 1/2 ) ~ T 2 ^(C^/2^)) = 4/// 2 *(V 3 (y) - T 2 ^(x)) 
= Q 3 - 3KQ - 2C 3/2 T n -x (C 3/2 /2v^) = 0, (3.42) 
where we define 

x = \{^f 12 = cosher), y = = cosh(^r), (3.43) 

and (C,Q) are the eigenvalues of the operators P ~ A and Q ~ — C T (see eq. ( I2.7ip ). 
One can easily see that this equation is an algebraic equation in ( and Q. Note that 
there are three branches, Q = Q m (() {m = 0, 1, 2), of the algebraic equation, 

F(t, Q) = {Q- Qo(C)) (Q - Qi(0) (Q - Q 2 (0) = o, (3.44) 

and this defines the three sheets of the Riemann surface. These branches are related by 
the following shift of the parameter r, 

^T^^r + 2mm (m = 0,1,2), (3.45) 

in the expression of eq. f)3.39p . since this shift does not change the coordinate Q. 
This algebraic equation has singularities (C*,Q*), 

= F(C*, Q*) = ^ = ^ , (3.46) 

at the origin (£*, Q„) = (0, 0) and at (£*, Q*) = (Cnj, Qn,j) of 

• / nn \ 2 / 3 _. / nir x 1 / 3 nn 

Q n ,j = u J n [2 cos — - - J , Q nJ = u J 2fA [2 cos - - - - J cos ^-' {3A7) 

with 

3 f 77,, 2q-l\n, neZ, (3.48) 

and j = 0, 1, 2. The geometrical meaning of singular points is that different branches of 
Riemann surface intersect at the points. Since the singular point at the origin does not 
change under the shift of r of eq. (I3.45p . this singular point joins the three sheets at the 
same time. On the other hand, the singular points given in (I3.47P only join two sheets 
at the same time. The reason is following: The invariance of the singular points (13.471) 
under the shift (13.451) is essentially given by the invariance of the factor cos ^ in Q n j, 
that is, 



cos — = cos 



n + (2q — l)m 



-7T 



(m = 0,l,2). (3.49) 
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Since the Z 3 symmetric p = 1 critical points satisfy 3 j g — 2 (=>• 3 { 2g — 1), the above 
condition is only satisfied by two different values of m(= 0, 1, 2). This means that the 
singular points only join two sheets at the same time. 
The branch points ((b,Qb) of the function Q((), 



appear when 



3R 



= F(C b , Qt 



2q — 1 | n 



dF(( b ,Q b 
dQ 



0^ 



dF(C b ,Q b 



jJCO J 



(.7 = 0,1,2) 



(3.50) 



(3.51) 



One can see that the singular points (I3.47P appear between the origin and the branch 
points in the Z3 symmetric manner (i.e. ( — > The simplest curve (the (1,1) cases), 



F(Q,C) = Q 3 -3/<Q-C a 



0. 



(3.52) 



is drawn in Fig. [2j The topology of the curves is also shown in Fig. [3j Note that there 
is a non-trivial handle in the algebraic curve. This comes from the cubic root in the 
curve. The appearance of non-trivial handle is explained in Appendix [0 where some 
basic properties of cubic roots are summarized. 




(H<0) 



Gix>) 



Figure 2: The algebraic curve of the 3-cut (1, 1) case and I = 2, which is given by F(Q,() — Q 3 — 
3^C ~ C 3 = 0. The branch points are at £ = 4 1 / 3 ^u/' (j = 0, 1, 2). The singularity is at the origin £ = 0. 



We also write the algebraic curve of the (2, 2) cases as 

^(C, Q) = (Q 3 - C 3 ) 2 + 3/ 8 CQ - 2f(Q 3 + C 3 ) 

with / = 3c/2 and 



(3.53) 



f^ 3 



3/ 



3/ 



2/3, 



A * 



3A 



3/, 



(3.54) 
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(a) 



(b) 



Figure 3: The topology of the curves are shown in the case of (a) (p, q) — (1, 1) and (b) (p, q) = (1, 3). 
The junction of three sheets is the singular point at the origin. The torus comes from the cubic-root cut. 
Two-sheet junctions in (b) are the singular points of (|3.47p . the number of which is nine. 

One can see that this solution also has three branch points. 

Next we consider the counterpart solution of "one-cut phase" in the two-cut matrix 
models, i.e. I = 0. We here adopt the expression (1C.12I) derived from direct analysis 
of the string equation, and leave the details in Appendix O The diagonalization of the 
matrix operators gives 

i~P~/izxO, -C T ~ Q ~v{z+ 1) x ft" 1 . (3.55) 

Here we put fi = t and z = K^d^ = ntd t (see eq. (13.21) ). They give three algebraic 
equations: 

F (n) (C,Q) =uo n Q -uj- n ( - n = (72 = 0,1,2). (3.56) 

Since the phase uj labels branches of the solution, the complete algebraic equation should 
be given by multiplication of these equations, 

F(C, Q) = F(°)(C, Q) ■ FW(C, Q) ■ F( 2 )(C, Q)=Q 3 - 3/iCQ - C 3 - ^ = 0. (3.57) 

This curve has no branch point but has three singular points which are at (C*,Q*) = 
(uji fi, —uj~^ )S) (j = 0,1,2). One can also see that these singular points only join two 
sheets at the same time. The curve is drawn in Fig. HI The topology of the curve is also 
shown in Fig. [5l 

Here in Fig. HI we draw a single connected cut which does not attach any branch 
points. This cut has a junction of three cuts. This kind of junction is always allowed in 
general k-th root (See Appendix [E]) . Since there is no branch point here, one can also 
push this cut away to infinity. 

Although this cut can be pushed away, this has a physical meaning as the eigenvalue 
condensation. There is a similar phenomenon in the two-cut matrix models [17,34]: In 
the one-cut phase of two-cut matrix models, the "two cuts" are completely connected 
with each other and form "a single cut" , which run along the whole real axes. Since the 
branch points disappear in this procedure, one can also push this cut away to infinity. In 
our three-cut case, since the cubic-root system admits a Z3 symmetric cut which has a 
three-cut junction without any branch points, it is natural to interpret that this system 
still have this kind of cut which indicates the eigenvalue density. 
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(H<0) (H>0) 

Figure 4: The algebraic curve of the 3-cut (1, 1) case and I = 0, which is given by F(Q,() = Q 3 — 
3/iCQ — C 3 — Z^ 3 = 0. The singularities are at £ — fJtuji (j = 0, 1, 2). There is no branch point, but we 
explicitly show the cut which indicates the eigenvalue condensation. This cut has a junction of three 
cuts, which is allowed in the cubic root. Since there is no branch point attached by this cut, one can 
also push this cut away to infinity. 




(a) (b) 



Figure 5: The topology of the curve is shown in the case of (p,q) = (1, 1). (a) is the case of finite \i 
and (b) is the case of \i = 0. 
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4 Conclusion and discussion 



In this paper, we have provided quantitative analysis of (p, q) critical points in the multi- 
cut two-matrix models based on an extension of the prescription developed by Daul, 
Kazakov and Kostov in the one-cut two-matrix models [27]. Right at the critical points, 
we identified the minimal construction of critical resolvents and potentials which can 
give (p, q) critical points of the multi-cut two-matrix models. The hermiticity of the Lax 
operators realized in the multi-cut critical points is identified. 

We have also studied off-critical behavior of macroscopic loop amplitudes. In this 
study, we identified an ansatz for non-polynomial parts of amplitudes in the Z& sym- 
metric background. This ansatz enables us to generate several formulas for off-critical 
macroscopic loop amplitudes: The intriguing one is the formula written with the Jacobi 
polynomials in the unitary cases p = q. 

Several future issues about our results are in order: 

1. In identification of the critical points, we only consider the system of equations 
(12.171) and (12.181) . In this investigation, there is no essential difference between the 
real potential models and the a/^-rotated models, in the sense of existence of Z& 
symmetry breaking critical points. As an independent check, one can perform direct 
evaluations using the Monte Carlo approach [49], since our real-potential solutions 
can avoid the difficulty of complex numbers in numerical analysis. 

2. Although we have found the formula for the unitary cases (13.361) . the expression 
in terms of algebraic equation is missing. Since underlying structure is multi- 
component KP hierarchy [33], they should be expressed by the algebraic equation 
of (kq,kp) order, F(x, Q) = [32,33]. 

3. From our analysis, the system of Douglas equation in the symmetric background 
admits one discrete parameter / in the expression (I3.36[) . The off-critical system 
is, however, controlled only by a dimension-full perturbation /i. This parameter / 
seems mysterious because one may expect that there are at most two meaningful 
geometries which respect the sign of /i@ The physical meaning of the parameter 
I, therefore, should be an interesting question to be investigated. 

4. Since our analysis only concentrates on the possible asymptotic (weak coupling) 
geometries of macroscopic loop amplitudes, the consideration about interpolation 
between different asymptotic geometries (which might be expressed by /) is nec- 
essary. This should be accomplished by non-perturbative analysis of the string 
equation as is in the two-cut one-matrix models [17,37]. 

5. Since we have known about macroscopic loop amplitudes in the Z& symmetric crit- 
ical points, it should be interesting to identify the Liouville continuum formulation 
corresponding to the Z& symmetric critical points. 

6. To see the relation to minimal fractional superstring theory, we need to calculate 
the amplitudes in the Z& symmetry breaking critical points. Since minimal frac- 
tional superstring theory is a natural generalization of bosonic and super string 

18 For example, see the two-cut case [17]. 
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theories, the formula could be probably written by the Jacobi polynomials of the 
ultraspherical sequence, Pn' a \z). 

7. One of the important features of the multi-cut matrix models is the charge of 
D-branes. In this sense, it should be important to see the annulus amplitude of 
macroscopic loop amplitudes [50,51]. To see the correlations among charges 
through annulus amplitudes of the matrix models is also an interesting problem to 
be studied. 
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A The relation between the resolvents and orthonor- 
mal polynomials 

An intuitive understanding of the relation between the resolvents and orthonormal poly- 
nomials is following: There is a useful expression for the general orthonormal polynomi- 
als [6]: 

a n {x) = -}j=(tet(x - X (n) )) = -i= ( e *M»-*(.>>\ 
AW .-i : ( d et( !/ -r ( „ ) )) w = -l : (e..Mv-, w ,) w . 

Here the expectation value is defined by path integral over the n x n truncated matrices 
X( n ) and F( n ) as 

(det(x-X (n) )\ = J- / d^ (n) dy (n) e- JVt " l '^Cn).y(»))(det(x-X (n) )) J (A.2) 

\ / {n) £(n) J \ / 

z (n) = n\ r n 

limit, 



with = n\ n«Lo hi- This implies the following powerful relation in the large N (~ n) 



A~ x , B~^(tr^), C~y, D-^tr^). (A.3) 

This relation has been argued and proved from several viewpoints [27] [29,52] [32,33]. 
With this relation, solving the eigenvalue problems (12.101) and (12.111) in the large limit 
is translated into obtaining the resolvents of the matrix model. 
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B Examples of critical potentials 

B.l Two-cut one-matrix models (p = 1) 

To compare our calculations with previous results in [14], we make a list of the critical 
potentials (of real coefficients) V^(x)@ Note that our (real-coefficient) critical potentials 
are related to the (complex-coefficient) critical potentials in [14] by replacing 6 — > ib and 
g — > i<7@ The parameter b corresponds to the critical point of eigenvalue (x* = 6), 
which one can always put zero, b = 0. The parameter g is an irrelevant perturbation and 
essentially does not change the system but is necessary to obtain the critical points of 
odd q (one can always choose any finite non-zero number). 

V x (x) = (-4b + 2b 3 -2g + b 2 g) x + - (-2 - 3b 2 - 2bg) x 2 + | x 3 + ^ (B.l) 

V 2 (x) = (2b -b 3 )x+ l - (-2 + 3b 2 ) x 2 - bx 3 + X —, (B.2) 

V 3 (x) = (-46 - 46 3 + 26 5 - 2g - 2b 2 g + b A g) x+ 

1 2 
+ - (-2 + 146 2 - 156 4 + 8bg - 8b 3 g) x 2 - - (b - 56 3 + g - 3b 2 g) x 3 + 

+ i(-l- 56 2 - 46^)x 4 + |x 5 + ^, (B.3) 

V 4 (x) = - (26 + 26 3 - 6 5 ) x + - (-2 - 66 2 + 56 4 ) x 2 + - (36 - 56 3 ) x 3 + 

+ ^(-l + 56 2 )x 4 -^ + ^, (B.4) 

y 5 ( x ) = I (-246 - 126 3 - 126 5 + 66 7 - 20g - 10b 2 g - 10b 4 g + 5b 6 g) x+ 
5 

+ — (_4 + 226 2 + 466 4 - 356 6 + 206^ + 406 3 # - 306 5 ^) x 2 + 
+ — (-86 - 646 3 + 846 5 - 10^ - 606 2 ^ + 756 4 ^) x 3 + 



15 

— (-2 + 366 2 - 1056 4 + 406^ - 1006 3 ^) x 4 
20 

— (-46 + 706 3 - lOg + 756 2 ^) x 5 + 



bq )x b + 1 . (B.5) 

15 10 V 7 40 V ; 

The critical resolvents are Qg defined by the operator £>( rcal ) (~ N~ l d/dx) as 

Q q (Z) = ^B^\*,Z) + B^ T (*,Z)y (B.6) 

19 One-matrix models can have a natural embedding in the two-matrix-model setup by putting V(y) = 
y 2 /2. The Gaussian integration over y for the two-matrix model partition function gives back to the one- 
matrix models. However note that the potential Vq{x) here is defined as the usual one-matrix potential 
Z = JdXe- NtrV( - x l 

20 Up to this replacement, we use the same notation as in [14]. The parameter b has a different meaning 
in the main text. 
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and 



Qi(z) = 



/(1-Z 2 )(1 + Z 2 



z 2 

(1-Z 2 ) 2 (1 + Z 2 

z 2 



(36 + g) 



(1 



^ 2 ) 2 (l+^). z \ 



<?4(Z) = - 



/ 

V 

(- (1 
V 

/ 



o 

(i-z 2 ) 2 (i + z 2 ) 
z~ 3 

Z 2 ) 3 (l + Z 2 ) 
2Z 4 

(i - zY(i + z 2 

2Z~ 5 



Z 2 

(1-Z 2 )(1 + Z 2 [ 

z 2 

(1- Z 2 ) 2 (l + Z 2 ) 



(3b + g)J 



z 3 



z- 1 

(56 + 2(7) 







/ 



(1 -Z 2 )\l + Z< 



■ z- 1 



2Z' J 

(1 - Z 2 f{\ + Z 2 ) 



■ z 







(1 -Z 2 ) 4 (l + Z 2 



Q 5 (Z) 



V 

/{l-Z 2 f{l + Z 



2Z 5 

2\5/ 



z- 



\ 



5Z 6 
(1 -Z 2 f{\ 



(76 + 5(7) 



• Z' 



2Z 4 

(1-Z 2 ) 4 (1 + Z 2 ) 
2^5 



(1-Z 2 ) 6 (1 



(56 



2^) y 



• Z 



\ 



hZ 7 

(1- Z 2 ) 5 (l + Z 2 ) 



5Z 7 5Z 6 

B.2 Two-cut two- matrix models (ft ^ 1) 

Here we turn on the coefficient 9 in the potential 

w(x, y) = Vi(x) + V 2 (y) + 9xy, 



(7b + 5g)J 



(B.7) 



(B.8) 



(B.9) 



(B.10) 



(B.11) 



(B.12) 



and put R* — 1. With suitable choices of the value of 9 or i?* (one of them is a free 
parameter), we can obtain critical potentials with smaller coefficients. For example, 
two-cut (4, 5) unitary critical potential with 9 = 1 is given as 



X 



-Ax 2 + 5670x 4 - 176400x 6 + 1157625x 8 (R, = 1/210), 



(B.13) 



and one can see that the coefficients are larger than those of the i?* = 1 case (1B.16|) . 
Also the critical point in eigenvalue space are chosen to be the origin (0,0). By 

making a shift of variable x — > ax + j3, we can eliminate the parameter 6 like in Appendix 

EU 

The unitary cases: 

The critical potential of the unitary cases V±(x) = V 2 (x) = V®\x) (i.e. (p,p + 1) unitary 
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minimal superstrings) are 



V®(x) 


= -2x 2 


4 

4 


= 6), 






V®(x) 


= 3x 2 - 


X 6 

3x 4 + -, 




(0 = 


40), 




V {A \x) 


= -Ax 2 


+ 27x 4 - 4x f 


X 8 

+ 8 ' 


(0 = 


210), 


V {5) (x) 


= 5x 2 — 


220x 4 + 70x 


6 - 5x s 


x w 
-\ 

+ 10 ' 


(9 = 1008), 


V {%) (x) 


= -6x 2 


6825x 4 


3250x 6 


555x 8 


_ 6x 10 + — , 
12 ' 


+ 4 " 


3 


+ 4 



(B.14) 
(B.15) 
(B.16) 
(B.17) 
4620). (B.18) 



The (p, q) cases: 

The potentials (V^\x), V^\y)) of general (PA) cases are 



(pA) 


= (2,3) 


(9 = 14) 


(PA) 


= (2,5) 


(9 = 39) 


(PA) 


= (2, 7) 


(9 = 76) 


(PA) 


= (3,4) 


(9 = 85) 


(PA) 


= (3,5) 


(9 = 154) 



+ £ + (B. 19 ) 

V} 5) (x) = 
Vi 2) (y) = 
V} 7) (x) = 



13x 4 


13x 6 


13x 8 


x 10 






8 


h Io' 


4 2 


\i 






2 


+ 4' 






95x 6 


57x 8 


57x 10 


19x 12 x 14 


6 „ 


2. 4 


5 


12 14 



V«» (I) = + 

2 vy; 2 4 g , 
Tr(5)/ . 55x 4 121x 6 llx 8 x 10 

2 vy; 2 26 
Z2 breaking cases: 

The Z 2 breaking critical potentials (g is a breaking parameter) are 

(p,q) = (2,3) (0=14) 

T ,f3i, v „ 15x 4 9ox 5 5x 6 ox 7 x 8 
^ J (x) = -3<?x + 3gx 3 + — 1 _ + -- + —, 

^ 2 (2) (y) = - 9 3 ) V + \ (-10 + Sg 2 ) y 2 - gy 3 + V - y 
(p,q) = (2,A) (9 = 25) 



(B.20) 



(B.21) 



(B.22) 



(B.23) 



(B.24) 
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Tr(4)/ , Agx 3 13x 4 _ 13x 6 12#x 7 13x 8 #x 9 x 10 

vfO/) = (<? - g 3 ) y + i (-13 + 3s 2 ) - ^ + v -. 

(B.25) 

Note that the critical potentials of cu 1 ^ 2 (= z)-rotated models are obtained by the analytic 
continuation, g — > ig. 

B.3 Multi-cut two-matrix models (Z& symmetric) 
B.3.1 Three-cut cases 
The unitary cases (p = q): 

The potentials of "unitary" cases V^x) = V 2 (x) = V®{x) are 

(B.26) 
(B.27) 
(B.28) 
(B.29) 
(B.30) 



V w (x) 


= - (1 - 4e - 5e 2 ) x 3 


ex 6 
+ 1P 


(0 = 


1 + 4e + 7e 2 , e ^ 1 


V^(x) 


8x 3 x 6 
= -3 + 6' <" 


= 12), 






V^{x) 


, 7x 6 x 9 

= 14x 1 

2 + 9 ' 


(0 = 


70), 




V^\x) 


200x 3 170x 6 


40x 9 


x 12 


(0 = 350), 


3 3 


9 + 


12 ' 


V {5 \x) 


910x 3 5005x 6 
= 3 " 6 + 


130x 9 - 


65x 12 
12 4 


™15 

■ — , (0 = 1638) 
15 V ; 



Note that, for the case of p = 1, we consider the potential of 7712 = p + 1 = 2. Then 
the critical operator of A is given as 

A(*, Z) = ^/r S 1 ~ Zk){ ^ ~ eZk) F(Z), (e = A^oo)). (B.31) 

From the construction, the vicinity of e = 1 is the p = 2 critical point. In this class of 
solutions, we need to be careful about critical potentials with odd degrees. In general, 
negative roots for the critical potentials V^ 2Z+1 )(:r) with I > 1, and the Fermi sea is 
filled from origin to some point around the first negative root. However for p = 1, the 
Fermi sea extends to negative infinity (See eq. (1B.26|) with e = 0). Hence we introduce a 



regularization parameter e by eq. (IB. 311) . This will make the matrix model well-defined 
and we will also adjust the normalization factor 9 accordingly. 

The (p, q) = (l,q) cases: 
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The potentials V\(x) = V x (x) and V 2 (y) of the (l,q) case are 

V} 2 \x) = I(-3 - 5e)x 3 + 

3 £ (B.32) 

I/ 2 ( y ) = -(l-9e)y 3 +^-, (0 = 3(1 + 3c), e^l), 

11 9 

K( 3 )( x ) = i e (i + I2e)x 3 + -(-5 - 8e)x 6 + — , 

3 6 g 9 (B.33) 

V 2 (y) = -(l-Ue)y 3 + e -^ rj (9 = 5 + 21e, e ^ 1), 
3 6 

y/ 4 >(x) = 1 (e + 15e 2 - 226 3 ) x 3 + +1 (7 + 23e + 33e 2 ) x 6 + i(-7 - lle)x 9 + ^, 
3 6 9 12 

V 2 (y) = ^(l-19e)y 3 + e -^, (9 = 7 + 38e, e ^ 1), 

(B.34) 

Since there is no (1,2) critical point in three-cut cases, the above (1,2) critical potential 
should result in the (1,3) critical point. 

The (p, q) cases: 

The potentials (V^ (x), V 2 ^\y)) of general (p, q) cases are 

r(3) 



(p, q) 


= (2,3) ( 


e = 26) 


(PA) 


= (2,5) ( 


9 = 69) 


(p, q) 


= (3,4) ( 


9 = 145) 


(PA) 


= (3,5) ( 


9 = 259) 



1 ^t 3 1 Q /tt«6 ™9 

2 vy; 3 g , 
Tr(5 ), N 161x 6 46x 9 23x 12 x 15 

3 9 12' 

6 9 

Tr(5)/ , 37x 3 407x 6 370x 9 37x 12 x 15 

(B.38) 



(B.35) 



(B.36) 



(B.37) 



B.3.2 Four-cut cases 
The unitary cases (p — q): 
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i(l-6e-7e 2 )x 4 + ^, (£ = 2(l + 3e + 5e 2 ), e^l). 



V^(x) 
V^(x) 
V^ 4 \x) 
V^\x) 



The potentials of "unitary" cases V x (x) = V 2 (x) = V®(x) are 
V w (x) 

4: 

.T 8 

18), 

100), 

(9 = 490) 



-3*< + ' T 
105x 4 15a; 



i 

-210x 4 
6435x 4 



8 x 12 

12' 



4 

175x 8 14x 12 x 16 
~~ 2 3" + To 7 ' 



1860x 8 + 



765x 



12 



45a; 16 x 20 
8 + 20 : 



4 4 
The (p, g) = (l,q) cases: 

The potentials V\(x) = V x (x) and V^iy) of the (l,q) case are 
V} 2 \x)= 1 -(-5-7e)x* + ^, 

(# = 5 + 13e, e^l), 



2268). 



(B.39) 
(B.40) 
(B.41) 
(B.42) 
(B.43) 



^(^ = 4(1-136)^ + ^-, 



i (3) ( 

V^\x) 
V 2 (y) 



\z/ 3) (x) = - (2 + 6e + lie 2 ) x 4 + -(-8 - lle)x 8 + 

I(l-20e)y 4 + ^, (0 = 2(4+15e), e^l), 



(e 2 + 25e 3 ) 



! ^ x 4 + - (22 + 61e + 60e 2 ) x 8 + —(-11 - 15e)x 
„ 8 12 



12 



(l-27e)y 4 + ^- 



ll + 54e, e^l), 



(B.44) 



(B.45) 



x 16 
~16' 



(B.46) 



Since there is no (1, 3) critical point in the four-cut cases, the above (1, 3) critical potential 
should result in the (1,4) critical point. 

The (p, q) cases: 

The potentials (V^\x), V^\y)) of general (p, q) cases are 



(p,g) = (2,3) 



38) 



(p,g) = (2,4) (0 = 65) 



(p,g) = (2,5) (5 = 99) 



v} 3) 


(x) = 




(y) = 


vl 4) 


(x) = 




(y) = 




(x) = 




(y) = 



19x 4 19x 8 x 12 

7 "§7 + "12 ' 

1% 4 , 



4 

■13x 4 + 



143x 8 13x 12 x 16 



13y 4 y s 



16 



33x 4 605x 8 99x 12 33x 



4 8 

33y 4 r 

'48 



16 



16 ^20 
20" 



(B.47) 



(B.48) 



(B.49) 
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287x 4 369x 8 41x 12 x 16 



205^ 4l/ 12 16 ( B ' 5 °) 

" ~~ 4 8~ + 12' 

455x 4 663x 8 377x 12 13x 16 x 20 



12 



(B.51) 



B.3.3 Six-cut cases 
The unitary cases (p = q): 

The potentials of "unitary" cases V x {x) = V 2 (x) = V®(x) are 

y( 1 )( x ) = I(l-10e-lle 2 )x 6 + ^, (# = 2(2 + 5e + 8e 2 )), (B.52) 
6 12 

V^\ x ) = -^ r + X —, (6 = 30), (B.53) 

V {3) (x) = 52x 6 - 4x 12 + — , (0 = 160), (B.54) 

18 

451 r 12 44r 18 r 24 

V^(x) = -748x 6 + - + (0 = 770). (B.55) 

The (p, g) = (l,q) cases: 

The potentials V\(x) = V 1 (x) and V^y) of the (1,5) case are 

Vf ) (x) = i(-9-lle)x 6 + ^, 

6 ^ (B.56) 

V 2 (y) = -(l-21e)j, 6 + ^-, (0 = 3(3 + 7e), 6^1), 

1 1 T 18 

^ (3) (x) = - (21 + 52e + 51e 2 ) x 6 + —(-14 - 17e)x 12 + — , 

y 12 12 18 

V 2 (y) = -(l-32e)y 6 + ^-, (0 = 2(7 + 24e), e^l), 
12 

v} 4) (x) = - (-19 - 71e - 141e 2 - 161e 3 ) x 6 + — (76 + 185e + 138e 2 ) x 12 + 

6 1 24 12 

+_(-19-236)x 18 + | I , (B.58) 

V^) = i(l-43e)/ + ^-, (0 = 19 + 86e, e ^ 1), 
o 12 

The (p, q) cases: 



(B.57) 



34 



The potentials (V}(x) : V^'iy)) of general (p, q) cases are 





= (2,3) 


(0 = 62) 


(p, q) 


= (2,4) 


{6 = 105) 


(P, q) 


= (2,5) 


{8 = 159) 


(p, q) 


= (3,4) 


(6 = 325) 


ip, q) 


= (3,5) 


(6 = 574) 



v{ 3) {x) 
Vi 2 \y) 



62a; 6 31a; 12 x 18 



V, 



(2 



(5 



V 



(2 



(4 



V 



(3 



V 



(3 



•J' 



2/) 



a; 



31y 6 / 2 
~6~ + "12 ' 

196a; 6 133a; 12 
+ 



18' 



7x 18 ^24 



3 v 12 4 

y 12 ' 



3 + 24 



;b.59) 



;b.6o) 



371x e 



30 



1007a; 12 265a; 18 53x 24 x 
f 6 24~ + 30 



y) = - 



x) 

y) 

x) 

y) 



3 fi I, 4 
53y 6 y 12 

"IT + "12 ' 

845a; 6 325a; 12 65a; 18 x 24 



:b.6i: 



IS 



1* 



24 



'B.62) 



3, J 

195^ 65y 12 y 

~~ 2 nr + is"' 

5863a; 6 5945a; 12 1927a; 18 41a; 24 x 30 



943y 6 Aly 1 ^ y 18 
~~ 6 6" + 18" 



12 



30 ' 
'B.63) 



B.4 Multi-cut two-matrix models (Z& breaking) 



We concentrate on the breaking critical points of the multi-cut matrix models which 
is characterized by the operators A and C T of 



— (1 - z k Y 
A(*,z) = ^r~ { - ] r(z), 



z 

c T (*,z) = v^ 1 1/ r(z)*- 1 



(l-Z^)4 

~ — r(Z), 



;b.64) 
;b.65) 



since they correspond to minimal fractional superstring theory [19]. The parameter g is 
the breaking parameter. The critical potentials of the a/^-rotated models are obtained 
by the analytic continuation, g (= Ci(*)) — > oj~ x g. We put R* = 1. 



B.4.1 Three-cut cases 

The unitary models (q = p + 1): 

The potentials (V^ +1 \x),V^\y)) of the unitary (p,p + 1) cases (unitary minimal frac- 
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tional superstring theory) are 



= (1,2) { 



V} 2 \x) = \e{-\ + 7e)gx 2 + -e(l + 12e)x 3 + -(-5 - 7e)gx* 

1 8 9 

, / r a % 6 , 9 X , 
+-(—5 — 8e )x ' 



6 



8 9 



V^iy) = +(-2 + 3e)gy + + 1(1 - 14e)y 3 - -egy* + 

2 3 4 6 



v} 3) (x) = 3gx 
Vi 2 \y) = gy + 



(p, g) - (2, 3) 13 ^3r 5gy4 y , 



3 i^i 



(# = 5 + 21e, e^l), 

8~ + 9 



13x 3 12gx 5 13x 6 gx 8 x 9 



.6 



4 + 6 ' 



26). 



(B.66) 



'B.67) 



The (p, q) = (l,q) cases: 

The critical potentials Vi(x) = V^\x) and V^d/) are 



V^ 3) (x) 



- (-e + lie 2 - 10e 3 ) gx 2 + - (e + 15e 2 - 22e 3 ) x 3 + 

1 1 

+- (7 + 19e + 25e 2 ) gx 5 + - (7 + 23e + 33e 2 ) x 6 + 

11 12 

+i(-7 - lOe)^ 8 + i(-7 - lle)x 9 + ^- + ^, 
= (-2 + 3e)gy + + 1(1 - 19e)y 3 - \egy* + ^, 

(fl = 7 + 38e, e^l), 



^ (5) (x) 



= --e (l - 19e + 51e 2 - 49e 3 + 16e 4 ) gx 2 + 
I 2 

+ - (e + 21e 2 - 96e 3 + 125e 4 - 51e 5 ) x 3 + e 3 (-l + 28e)#x 5 + 



(B.68) 



+- (e 3 + 34e 4 ) x 6 - - (22 + 93e + 188e 2 + 192e 3 ) #x 8 + 

+- (-22 - 105e - 234e 2 - 255e 3 ) x 9 + — (33 + 95e + 88e 2 ) gx n + 

1 1 
+— (33 + 103e + 102e 2 ) x 12 + — (-11 - 16e) gx 1 

15 V ; 17 18 



,14 



V2{y) 



-2 + 3e)gy + \eg 2 y 2 + 1(1 - 29e)y 3 - \egy" + ^, 

(0 = ll + 87e, e^l), 



(B.69) 



B.4.2 Four-cut cases 

The unitary models (q = p + 1): 

The potentials (V^ +1 \x),V^\y)) of the unitary (p,p+ 1) cases (unitary minimal frac- 
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tional superstring theory) are 
(p,g) = (l,2) 

Vi (2) (x) = - (4 + 8e + 15e 2 ) #x 3 + - (2 + 6e + lie 2 ) x 4 - -(4 + 5e)#x 7 

O — ( 

1 . . Q fifX 11 

+ 8 ( - 8 - 11 ^ 7 + TT + l T2' 

^ 2 (1) (2/) = ("3 + ^)gy + j^V + ^1 - 20e)y 4 - -e^ 5 + |y 8 , 

(0 = 2(4+15e), e^l), (B.70) 

(p,g) = (2,3) 

if )(x) = -^x 3 - 13x 4 + ™Zj + ^ - ^ - ^ + ^ + ^ 
1 K J 3 7 8 11 6 15 16' 

(0 = 65). (B.71) 

The (p, g) = (l,g) cases: 

The critical potentials Vi(x) = V^\x) and V^y) are 

V} 3 \x) = - (3e 2 - 28e 3 ) gx 3 - - (e 2 + 25e 3 ) x 4 + - (22 + 54e + 49e 2 ) gx 7 + 

3 1 2 1 7 1 

+ - (22 + 61e + 60e 2 ) x 8 + —(-11 - 14e)#x n + —(-11 - 15e)x 12 + 
8 11 12 

^gx { x 



15 16' g 
V 2 (y) = (-3 + Ae)gy + -e^ 2 + 1(1 - 27e)y 4 - 7 -egy 5 + 

(0 = ll + 54e, e^l), (B.72) 

V^x) = l e 2 (3 - 40e + 45e 2 ) gx 3 + -e 2 (-2 - 60e + 95e 2 ) x 4 

\ 1 
— (14 + 51e + 91e 2 + 84e 3 ) gx 7 + - (-28 - 117e - 232e 2 - 228e 3 ) x 8 

7 1 1 8 2 

(49 + 124e + 99e 2 ) gx 11 + — (49 + 134e + 114e 2 ) x 12 (7 + 9e)gx 15 



11 

+ ±(_14_l 9e)x i6 , 9X 



W 15 



16 v ' 19 20' g 

V 2 (y) = (-3 + 4e)gy + ~e<?y + i(l - 34e)y 4 - ^e^ 5 + 

(0 = 14 + 85e, e^l), (B.73) 

B.4.3 Six-cut cases 

The unitary models (q = p + 1): 

The potentials (V^ (x),V 2 (y)) of the unitary (p,p + 1) cases (unitary minimal frac- 
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tional superstring theory) are 
(p,g) = (l,2) 

' V} 2) (x) = - (21 + 44e + 40e 2 ) + - (21 + 52e + 51e 2 ) a; 6 - —(7 + 8e)gx u + 

1 . . in OX 17 X 18 

+—(-14 - 17e )x 12 + + — , 

12 v ; 17 18' 12 

^2%) = (-5 + 6e)gy + ^e«?y + i(l - 32e)y 6 - ^ 7 + 

2 b 7 Iz 

(# = 2(7 + 24e), e^l), (B.74) 

(p,g) = (2,3) 

y (3) f ^ = _287^ 5 _ 196^ 36% u 133*^ _ 41p 17 _ 7* 18 g 23 x^_ 
1 K ' 5 3 11 4 17 3 23 24' 

T/(% \ , n 9 2 7 6 u 9 7 . 2/ 

^ 2 u/) = #y + -g-j/ - V - —j/' + — . 

(0 = 105). (B.75) 

The (p,q) = (1,<?) cases: 

The critical potentials Vi(x) = V®(x) and V 2 {y) are 

V (3) ( x ) = - (-19 - 58e - 100e 2 - 110e 3 ) gx 5 + - (-19 - 71e - 141e 2 - 161e 3 ) x 6 + 
5^ ^6 

+ — (76 + 172e + 121e 2 ) gx 11 + — (76 + 185e + 138e 2 ) x 1 

1 1 or 23 r 24 

+^(-19 - 22e)^ 17 + ^(-19 - 23e)x™ + ^- + * 

V 2 (y) = (-5 + 6e)gy + ^eg 2 y 2 + ~(1 - 43e)y 6 - ^egy 7 + 6 



,12 



2 * ° 6 V 7 12 

(#=19 + 86e, e^l), (B.76) 

y/ 4) (a;) = - (6 + 24e + 60e 2 + 120e 3 + 245e 4 ) gx 5 + 
5 i 

+ - (3 + 15e + 45e 2 + 105e 3 + 203e 4 ) x 6 - 
3 

-— (56 + 189e + 264e 2 + 154e 3 ) gx n + 
+ — (-224 - 819e - 1224e 2 - 754e 3 ) x 12 + 

T i 

+ Yj ( 78 + 180e + 119e2 ) g^ 7 + - (52 + 126e + 87e 2 ) x 18 - 

_±( 6 + 7e ) gx ™ + _L(_ 24 - 29e)x 24 + 9 — + — , 
23 V Jy 24 v ; 29 30 ' 12 

7 2 (y) = (-5 + 6e)gy + ^eg 2 y 2 + ~(1 - 54e)y 6 - ^egy 7 + 

(# = 3(8 + 45e), e^l), (B.77) 
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C Direct evaluation of three-cut string equations 



In this appendix, we make a detailed analysis of the string equation in the three-cut matrix 
modelgj^J based on the formalism of fc-component KP hierarchy [53] [33] . The motivation 
behind this analysis is to check that our ansatz (13. 23ft corresponds to a special solution of 
Lax pairs. On the other hand, the complete analysis also provides useful hints for other 
solutions in the multi-cut matrix models. 

The realization of symmetric (p, q) critical points have been shown in eq's (12.711) 
and (12.731) . Here we drop some irrelevant factors and the following simplified notation is 
used: 

P = + #i d?- 1 + ■ ■ ■ + Hp, Q = r- 1 ^ + H x d«~ l + --- + H 9 . (C.l) 

The coefficient matrices Hi and Hi are all real functions. Then the Lax operators T and 
L are defined as P = (Tlf) with [P, T] = [P, L] = and T k = 4, and 

oo oo 

r = r + s n (t) d~ n , L = I k d + J2 x n(t) d~ n . (C.2) 

n=l n=l 

In the following discussion, we put k — 3 for simplicity. The relation T 3 = J 3 can be 
easily solved if derivative terms are neglected: 

£(diag) = _r j2 {S ni ,S n2 ,S n3 } + (derivative terms), (C.3) 

ni+n2+ri3=n 

with So = T. The bracket { , , } is a totally symmetric product which is defined by 

oo g 

= E i S n,S, 2 ,S l3 } (C.4) 

Here we define 

F (dia g ) ^ Ijr r ^^ F (os) = I[r 2 , [t,f]], f = F (diag) + F (off) . (c.5) 

3 3 

In the KP hierarchy basis (r — > Q), they are truly a diagonal part and an off-diagonal 
part of the matrix F. 



C.l The p = 1 cases 

Since we impose the symmetry here, the KP operator P corresponding to the matrix 
operator A is given in eq's (12. 711) and (12. 73ft as 

P = Yd + H = (TL), (C.6) 

with 

H=\ h\, fx + h + h = 0- (C7) 



21 Thc calculations in this appendix have been carried out mostly with the help of Mathematica™ . 
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Then the equation [P, T] =0 gives 



^ off ) = -- [r 2 , [H, S'n-i]] + (derivative terms). (C.J 



The operator Q for the (l,q) case can be written [48,54] as 



Q = Y j nt n (T~ 1 L n ~ 1 ) + , (C.9) 

n=l 



with (q + = 1. Then the string equation can be easily derived 22 ! from 



9+1 



= ^nt„[P,(r" 1 J L"- 1 ) + ]. (CIO) 



n=i 



One can easily see that q = 2 (mod. 3) is not allowed as critical points. For example, in 
the first simplest case (1, 1), the string equation is written as 

= T0HT+ [H 2 ,T] +t[H,T 2 ] (C.ll) 



with t\ = t. Below, by imposing the information of the matrix model (1C.7|) . we solve this 
string equation at the first order of string coupling. 

C.l.l The case of (1, 1) 

There are four asymptotic solutions: The first one is generalization of "one-cut solution" 
in two-cut matrix models: 

P = Td, Q = r- 1 (9 + t), (f 1 = f 2 = 0). (C.12) 

With a shift of z, this corresponds to our solution (13.361) with the choice of / = 0. The 
others (i = 1,2,3) are essentially expressed as 



with 



Pi = Td-tMi, Qi = T- l d -tMi, (C.13) 



M x =| -2 j , M x = (-2 | , (C.14) 

M 2 = ( 1 j , M 2 =^l | , (C.15) 

M 3 = | 1 j , M 3 = (-2 | . (C.16) 



J See Appendix D in [33], for example. 



40 



They all satisfy Mf = Mj. Note that the forms of Pj and Q i (i = 1,2,3) are different 
from each other but eigenvalues are the same. With a shift of z, this is equivalent to our 
ansatz (13. 23j) and the solution (13.361) of / = 2. 

We also note that the Douglas equation, [P, Q] = g s t r h, requires the following first 
order correction of string coupling: 

P = Td - tMi + g str jJi + 0(g 2 str ), (C.17) 

with 




(C.18) 



This kind of first order correction cannot be realized by the ordering problem of z and t 
like in the two-cut cases ( I3.17p . 

C.1.2 The case of (1,3) 

We consider the background {b\ = t,b 2 = t 2 , b 3 = 0, 6 4 = 1/4). There are five asymptotic 
solutions: The first one is a generalization of "one-cut solution" in two-cut matrix models: 

Po = Td, Q = r" 1 (9 3 + t 2 9 + t), (/ 1 = / 2 = 0). (C.19) 

Other three solutions (i — 1, 2, 3) are essentially expressed as 

Pi = Yd + fMi, Qi = (d 2 + fd-f + t 2 ) (T^d + fMi) , (C.20) 

with 

3t = (h - \f)f- (C.21) 

These solutions precisely correspond to the ansatz (13.231) . 

The last solution is different from the ansatz (I3.23p . and can be expressed as 

P 4 = Td + H, Q A =(Td + H) 2 d = P\d. (C.22) 
The functions f\ and fa in H (see (IC.7D ) are related as 

t = \hh{h + h), h = l -ul + hfa + f 2 2 ). (C.23) 

It would be of interest to make further study of these solutions. 
C.2 The case of (2, 2): p = 2 

The KP operator P which corresponds to the matrix operator A is given in eq's (12.711) 
and fT2T73]) as 

P = Td 2 + H 1 (t)d + H 2 {t) = (TL 2 ), (C.24) 
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with 



fi 



9i 



#1 



92 



fi + / 2 + fs = 0. 



(C.25) 



,93 



Then the equation [P, T] =0 gives 
1 



c(off) 



[r 2 , ([H h S n ^} + [H 2 , S n ^ 2 ]\] + (derivative terms). (C.26) 



The operator Q and the string equation for the (2, 2) case can be written as 

(T~ l L 2 )+, \YL\T- l L 2 ] 



Q 



9strl-i 



(C.27) 



There are also several equivalent solutions like the 2nd, 3rd and 4th solutions of the (1,1) 
and (1, 3) cases. Here we only list representative solutions. 

The first solution is the counterpart of the one-cut solution: There are three equivalent 
expressions. One of them is given as 



/ 



t 2 

2. 



\ 



d- J - 



t 2 
2. 



d 



f 







Q 



( 



\ 



o 



t 2 

2. 







d- J - 



f 



A 2 \ 



/ 



(C.28) 



The parameter / is fixed by the Douglas equation [P, Q] = g s t r h- This case is also 
related to our solution (13.361) of / = 0. 

The second solution corresponds to our ansatz (13.231) : There are also three expressions 
and one of them is given as 



/ 



-J 



( 



d-lf 



d+-f 

r 



V 



d--f 



(C.29) 



/ 



The parameter / is fixed by the Douglas equation [P, Q] = g s t r Iz- This case is also 
related to our solution (13.361) of I = 2. 
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D The Jacobi polynomials 



Here we write some basic facts about the Jacobi polynomials, which we have used in this 
paper. Some standard reference is [55], for example. 

D.l A basic definition and the differential equation 

The Jacobi polynomial Pn a,l3 \z) is a solution of the hypergeometric equation, 

\l-z 2 )^ + (p-<*-(a + P + 2) z)A + n(n + a + + i^Pj*® (z) = 0, (D.l) 
and is defined by 

l 1 'in ^yn r -. 

1 - z) a+n {l + zf +n . (D.2) 



2 n n\ dz n 
They are orthogonal polynomials with respect to the following inner product: 

1 dzPM(z)pM(z)(l - zT(l + zf = 2 " P+1 nn + a + inn + P + 1) 

azi> n \Z)f m [Z){L Z) [L + Z) 2n + a + p + 1 nW{n + a + p+1) <W 

(D.3) 

D.2 Connection with other polynomials 

The Jacobi polynomials are related to other famous polynomials. The Gegenbauer poly- 
nomials, Cn\z), and the Legendre polynomials, P n (z), are 

p(A-l/2,A-l/2) (z) = (A+l/2) w C £ ){zl p ^ ){z) = pn{zy (D 4) 

The Chebyshev polynomials of the first kind, T n (z), and the second kind, U n (z), are 

p(-l/2,-l/2)/ n _ ( 2n ) ! rp / \ p(l/2,l/2)/ N (2ra+l)! 

n { } ~ 2 2n (n\) 2 { h n { 1 ~ 2 2n+1 [(n + l)!] 2 ^ >' { 1 

They are also called the ultraspherical polynomials P^ a,a \z). The twisted type of ul- 
traspherical polynomials, p( a '~ a \z) [56], can include the Chebyshev polynomials of the 
third kind, V n (z), and the forth kind, W n (z): 

P »"' A1,2,(2) = P »" 2 '" /2,(Z) = W " {Z) - (D - 6) 

So we can see that W n (z) = (— l) n V n (— z). 

D.3 Useful relations 

Some useful relations are in order: 

Reflection relation: 

Pi a ^(z) = (-irP^\-z). (D.7) 
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The normalization: 

P (^) (1) 

The leading coefficient: 



n + a 
n 



P(^)(-1) = (-1)' 



n + (3 
n 



T(2n + a + p + 1) 
2 n n\T(n + a + (3 + 1) 



z n + - 



Derivative formula: 



d pM(~\ - n + ® + P+ 1 p (a+l,P+l) ( x 



D.4 Proof of the solution f T3736T > 

Here we check that our formula (I3.36P 

■2l-k 2l-k 



k-l 



I 21- k 21 — k\ I ii 



(z-iy- l (z+iy, 



generally satisfies eq. (13. 6p 



?*$(*K«(*0 -p^(z)7rp(z) = const., 

in the case of p = q. For sake of simplicity, we neglect normalization factors, 
need to show is that the derivative of eq. (ID. 121) vanishes, 

*}(z)Ze(z)-Q{z)irp(z) = 0. 
We first calculate n'Uz) £,p(z) as follows: 



,2l-k 2i-fci 



f)2 pV k ' fc 
u z r p-\ 



(,)x(/(,-l) l (z + l)^ + 



[ 21— k 21— k \ 

+ 2d z Pl_1 ' k \z) xd z 



"z-l) l (z+l) 



k-l 



I 2l-k 2l-k \ 
fc ' k't 



fP^X' ^\z)xd 2 z (y(z-l) l (z+l) k 

( 2l-k 21-— k \ 



k-l 



I 2l-k 2l-k \ 

pL k ' k \ z 



+ 2[z+- 



l k — l\„ _2i^fc) 



k k 



d P k 



I (k-l 

k\ k ) D ( i r ± 



/_ k 



L( k-l 
k\ k 

z 2 - 1 



+ p(p- 1) 



/ 2l-k 2l-k \ / 2l-k 2l-k \ 
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Here we use the following relations, 

{/(z - lf-'iz + 1)'x(/(j- \)\z + I)'"' = : 2 - 1, (D.15) 



J (k-l 



and the differential equation (ID. II) , 



(D.18) 



Since the calculation of £,p(z)iip(z) is similar to that of n'l(z) £$(z) by exchanging the two 
indices: 

s " — ■ (D - 19) 

and one concludes 

tf(*0 <(*)&(*)■ (D.20) 



This is the equation (1D.13|) . One can calculate the constant in eq. (1D.12j) and fix the 
normalization of eq. (I3.36|) by using eq's (1D.8j) and (1D.9h . 



E Topology of cubic roots 

Here we summarize basic properties of cubic roots of analytic functions and the topology 
of associated Riemann surface, which serves as the simplest prototype of general k-th 
root cases. In the case of cubic root of a single complex variable, 

v 7 ? -> f/(e 2m zy = uj 1 \/7 1 = 0,1,2,---), (E.l) 

we can define the charge q of the branch point 2 = 0, according to the power of the 
fundamental phase I, as q = I mod 3. One can have branch cuts in the complex plane 
which connect branch points with vanishing total charges, q! + q 2 + • • • = mod 3. For 
cubic roots of an analytic function, there are two types of branch cuts: 

(i) Type I branch cut connects two branch points of opposite charges. One concrete 
example is given by 

h{z) = y{z-a){z-b)*. (E.2) 

Here the branch point z = a carries charge q = 1, and the branch point z = b carries 
charge q = 2. The geometry of the associated Riemann surface is given in Fig. [6(a). By 
cutting the three complex planes along the branch cuts and gluing the six edges according 
to the identification rule, we can see clearly that the topology of the type I branch cut is 
given by a simple junction of the three sheets, as shown in Fig. [6(b). 
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Figure 6: The first type cubic cut of finite length, a) How to connect the three sheets, I, II and III. 
The symbols a, /?, 7 indicate the way of connection, b) After cutting the cuts out and gluing them, one 
can observe that the topological role of the first type cut, which is a junction of the three sheets. 



(ii) Type II branch cut connects three branch points of identical charges. One concrete 
example is given by 

f 2 (z) = ^/{z-a){z-au){z-au 2 ) = v^ 3 - a 3 (E.3) 

Here all three branch points ( = a,au,au 2 carry charge q = 1. The geometry of the 
associated Riemann surface is given in Fig. [7(a). By cutting the three complex planes 
along the branch cuts and flipping the resulting complex plane into a hexagon (including 
the point at infinity), one can transform the identification rule in Fig. [7(a) for three 
complex planes into a gluing rule for three hexagons. Hence, the topology of the type II 
branch cut is the same as that of a torus (Fig. 0(b)). 




Figure 7: The second type cubic cut of finite length, a) How to connect the three sheets, I, II and III. 
The symbols o;,,,9,-,7, indicate the way of connection around the branch point i (= 1,2,3). b) Cut the 
cuts out and open them up. After gluing the sheets, one can observe that the topology is torus. That 
is, the topological role of the second type cut is a junction of the three sheets with a torus. 
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